Local density and the RPA corrections in 
charge current quasielastic neutrino on 
Oxygen, Argon and Iron scattering 

Krzysztof M. Graczyk* 

Institute of Theoretical Physics, University of Wroclaw 
pi. M. Borna 9, 50 - 204 Wroclaw 
Poland 

February 8, 2008 



PACS: 25.30.Pt 



Abstract 



Numerical computations of cross sections for quasielastic charge current 
scattering of neutrino on Oxygen, Argon and Iron in Local Density Approx- 
imation (LDA) are presented. We consider processes for a few GeV neutrino 
energy. We include corrections from nucleon re-interaction in nucleus de- 
scribed by relativistic Random Phase Approximation (RPA). We adopt the 
relativistic Fermi gas model of nucleus with and without taking into account 
the effective mass of nucleons. 



1 Introduction 

The study of oscillations phenomena of neutrinos became one of the 
most interesting topics of particle physics. Experiments such as K2K, 
MINOS, MiniBoone, ICARUS collect or will collect data which allow 
to establish parameters of oscillations. In the analysis of experimental 
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results a good knowledge of theoretical predictions of neutrino-nucleus 
cross sections is crucial. In long baseline and atmospheric neutrino 
experiments one of the most important become analysis of the data 
for a few GeV neutrinospQ. In this energy domain the quasielastic 
contribution is important. 

The calculations of cross sections require a study of influence of 
nuclear effects. The scattering amplitude can be combined from two 
parts^. One is described by weak interaction vertex of scattering neu- 
trino on free nucleon and the other which is evaluated in the framework 
of many body theory and describes the model of nucleus. The weak 
interaction vertex is described very accurately by Fermi theory a few 
GeV neutrino. A consideration of the model of nucleus described by 
Fermi gas (FG) of nucleonsjS] is one of the most popular and simple 
starting point in discussions of nuclear effects in that energy regime. 
In the Fermi gas model momenta of nucleons are uniformly distributed 
inside the Fermi sphere whose radius is called Fermi momentum - hp. 
It is connected with baryon density inside nucleus. For the first time 
this simple approach was applied in the scattering of electrons on nu- 
clei pQ-[Z|. Experimental results obtained from these processes allowed 
to establish average Fermi momenta for given nuclear targets. Suc- 
cess of this approach caused for application of it in neutrino-nucleus 
interactions [HI El EB G3] • 

As was mentioned above the Fermi gas model is the good starting 
point for further nuclear effects discussions. One can evaluate this by 
consideration of some additional effects such as binding energy, effec- 
tive mass of nucleons, final state interaction and many others. However 
the most important seems to be taking into account the collective be- 
havior of nucleons inside the nucleus. It can be done in several ways 
e.g. by using static potential, which gets the quite easiest way to apply 
it or by application of the methodology of quantum many body field 
theory. 

One of the most consistent ways of describing of the collective prop- 
erties of nucleons inside nucleus is the relativistic Mean Field Theory 
(MFT)p2]. The ground state of nucleus is described by relativistic 
Fermi gas of nucleons which has effective mass introduced by self- 
consistency equation for given target. The strong interaction between 
nucleons is introduced into the model by consideration of virtual par- 
ticles which are exchanged between nucleons. They are described by 
the effective interaction lagrangian which contain of vector and scalar 
fields. The MFT allows to solve this theory. In the electron-nucleus 



2 



scattering, neutral a and vector u are considered ^H]- The contribution 
from a field is fully manifested in the effective mass (the appearing of 
this field lead to the above mentioned self-consistency equation). In 
the case of the charge current neutrino-nucleus interaction, pions and 
mesons rho give contribution to scattering amplitnde|14| I15j. Correc- 
tions described by effective lagrangian are calculated by using standard 
tools of quantum perturbation field theory. It gives rise to infinite 
number of diagrams. It is impossible to take into account contribu- 
tions from all diagrams. One may chose special class of them which 
correspond to one particle - one hole excitations. This approach is 
called relativistic Random Phase Approximation (RPA) [16J. 

Experimental measurements of the electron-nucleus scattering show 
that charge distribution in nucleus is not uniform[TTjlT£]. which means 
that baryon density is not uniform either. The charge distribution can 
be calculated in the framework of the relativistic MFT by solving the 
Dirac-Hartree self-consistency equations^H]- Solution of it is consis- 
tent with experimental predictions. It allows us to use the experimen- 
tal results and apply them to the formalism under consideration. 

The measurements of the scattering electron on nucleus allowed to 
fit charge density distributions inside nucleus. The fitting is done by 
assumption that the simple models of nucleus such oscillator model 
in the case of light nuclei e.g. Oxygen and Fermi model in the case 
of heavier nuclei e.g. Argon, Iron. The inclusion into the model of 
nucleus of the local profile of baryon density (Local Density Approx- 
imation - LDA) is the treatment which makes the model of nucleus 
more realistic. The effect of the local density were taken into account 
not only in discussions of electron nucleus data and in neutrino matter 
interactions[20J-[23j. 

In this paper we want to develop consideration of previous paper [24|, 
where we adopted relativistic MFT. We presented the analytical solu- 
tion of the RPA Dayson equation and discussed the RPA corrections 
in scattering of neutrino on nucleus. Here we are going to extend our 
consideration by application of Local Density Approach (LDA) in the 
cases of these few targets. We will compare results for the case of LDA 
with calculations which were done by assumption the constant baryon 
density. The second case is described by average values of Fermi mo- 
menta and effective masses for given targets. They were calculated 
from experimental charge density distributions[18j. 

We choose nuclei sO 16 , is^ 40 and which characterize tar- 

gets of neutrino experiments: Super Kamionkande, Icarus and Minos. 
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As mentioned above, we are going to compare differential and total 
cross sections for the nuclei with constant and local density. The com- 
parison be done in various levels of the model (Fermi gas with or 
without RPA and with or without effective mass). 

The paper is organized as follows. A short description of the for- 
malism is presented in section 2. In section 3 numerical computations 
of cross sections are presented and the results are discussed. In this 
section we also summarized. 



2 The formalism 

We consider the process: 

+ N(A, Z) -> n~+p + N(A - 1, Z) 

The cross section per one nucleon in the laboratory frame for the above 
scattering is the following: 

d*o _ G|cos^ c |,1 Im(Vn ^ (1) 



d\q\dqo 16tt 2 pfE 2 

where pp = k F /37T 2 , q = k — k' is the transfer of four-momentum, k, 
k' are neutrino and lepton four momenta. L^ v is a leptonic tensor: 



8 (yk^k v + k ^k v fffiukak i i£^af3k k^ 



The ± sign depends on the considered process (neutrino/antineutrino). 
Nuclear properties of nucleus are included in the polarization tensor 
U^ u . This tensor is evaluated by the QFT techniques ^H] and can be 
written as the sum of two contributions. One of them is the free Fermi 
gas and second the RPA correction: 

n^(g) = n^ ee (g) + n^ PA ( g ), 



HSU*) = -i/ tSi^ ( a fg(p + jTOWH)) • (2) 



d 4 p 

By elementary weak charge current nucleon-nucleon is described 
and expressed by means of the form factors |25| : 

\a av n 

r a (g») = iM^h Q + ^(gj)-^ + ( 3 ) 
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We omit G p {q^) term since its contribution to v e and cross sections 
at E v ~ lGeV is negligible. The value of parameters which character- 
ize the form factors are presented in the appendix. 

By the A FG (p) in © we denote the propagators of nucleon in the 
free Fermi sea: 

A FG ( P ) = (jf + M*) ( ^_^ 2 + i£ + |^*(po - - p)) . 

(4) 

M* is the effective mass, calculated from the self-consistency equation 
for given Fermi momentum k F |12j : 

#f = (4 + M* 2 ) 1/2 . 

In our calculations we will separate the case with the free nucleon mass 
M* = M = 939 MeV and with the effective mass. 

The RPA part of the polarization tensor is calculated by consider- 
ation of contributions from one particle - one hole excitations. It leads 
to the correspond Dyson equation whose solution gives the n^p^(for 
more details see |24|). 



2.1 Local Density Approximation 

We are going to compare cross sections for the case of the model of 
nucleus with constant baryon density with the case of the model of nu- 
cleus with local baryon density. The constant value of baryon density 
for given nucleus is evaluated by calculation of average value of Fermi 
momentum for given density distribution 

< kF > = i^iMMO , (6) 

k Fp (r) = pr 2 |p(r), (7) 



k Fn (r) = pr^(r). (i 
A = J d 3 rp(r) is a atomic number, p(r) - charge density profile[18j. 
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The results for the local baryon density are obtained by the in- 
tegration of the cross section per one nucleon with weight given by 
density. 

r cutoff 

dviocal density = 47T / drr 2 p ntP (r)da(k F (r) , M* (r)) | pcr nucleoli , 

J 

(9) 

where the index of n or p in p corresponds to neutron or proton density 
distribution. 



Nucleus 


<k F > [MeV] 


< M* > [MeV] 


Oxygen 


199.21 


690.78 


Argon 


216.88 


631.37 


Iron 


216.91 


634.84 



Table 1 ! The average values of Fermi momenta and corresponding effective masses calculated from 
charge distributions for nuclei. 

In presented approach, for each values of integrating variable r the 
Fermi momenta fcjr(r) and effective masses M*(r) are calculated. It 
is done numerically. The reconstructed distributions of the effective 
mass inside nuclei are presented in fig. ^ We also calculate the average 
Fermi momenta for Oxygen, Argon and Iron as well as effective masses. 
The results are presented in tabled 



3 Numerical results and discussion 

We show a set of plots which present the differential cross sections for 
1 GeV neutrino energy. We calculated also the ratio of total cross for 
the model of nucleus with constant density to cross sections calculate 
in LDA at 1 GeV neutrino energy (tableEJ). In the end we computed 
the total cross sections for the three investigated targets. It was made 
by taking into account all considered nuclear effects (LDA, RPA cor- 
rections and the effective mass). 

In the table [21 we made separation of several cases. We computed 
ratio of cross sections for the free Fermi gas, the Fermi gas with effective 
mass, the Fermi gas with RPA without and with the effective mass. 
In the case of the model of nucleus with constant density the average 
Fermi momenta and the effective masses from table ^ were applied. 
As can be seen the differences between approaches are in order of 
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Nucleus 


FG 


FG (M*) 


FG + RPA 


FG + RPA (M*) 


Oxygen 


0.997 


0.994 


0,997 


0.993 


Argon 


0.994 


0.986 


0,994 


0.989 


Iron 


0.994 


0.987 


0,995 


0.990 



Table 2: The ratio of total cross sections calculated for nucleus with local density to cross section for 
the nucleus with constant density. The calculations were done for 1 GeV neutrino energy. FG means the 
free Fermi gas model. The appearing of M* means of application of the effective mass. 



1.5%. It does not depend on how sophisticated complications the 
model of nucleus is. One can also notice that the differences are slightly 
bigger for argon and iron. These two targets are heavier and charge 
distributions for them were extracted in the different way. 

Make our discussion more precise we present details comparison of 
differential cross sections. In the figure El and 01 the differential cross 
sections for scattering of 1 GeV neutrino on Oxygen, Argon and Iron 
are presented. It was calculated for the Fermi gas model of nucleus 
with and without adoption of the effective mass. The differences be- 
tween approaches with constant and local baryon density of nucleus 
are minor. One can notice that only the maximum of the curve in the 
LDA is little cut. The application of effective mass lowers the pick of 
cross sections by about 30 % but also does not change the differentiae 
between results for LDA and constant baryon density. 

In last three figures we consider the Fermi gas with corrections 
from random phase approximation. RPA corrections are calculated by 
counting infinite number of one loop digrams^!]. We consider charge 
current reactions with mesons p and pions contributions into RPA 
polarization tensor. We applied standard values of coupling constans 
for these fields. The short range correlation between nucleons is applied 
by using Landau-Migdal parameter g' |14j . In our calculations we 
adopt g' = 0.7. 

In figure |3] we present calculations with the RPA corrections and 
M* = M = 939 MeV. Similarly, as was noticed in the previous case 
(figEJ, only small differences in the maximum are observed. The in- 
clusion of the effective mass in the RPA calculations dose not change 
the LDA effect (fig. When we compare the figure 01 with fig. 0]and 
figure El with figre[£]one can see that the RPA effect decrease the cross 
sections, especially for small transfers of the energy. 

We conclude that local density effect could be replaced by calcu- 
lation of equivalent average parameters of the model which make the 
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numerical calculations faster. The application of the approach does 
not depend on complication introduced dynamical model of reaction. 



A Appendix 

• Form Factors 

In our computations we used dipol form factors PH] with: M A 
1.03 GeV , My = 0.71 GeV 2 , fi = 4.71 and g A = G A (0) = -1-26. 



Charge density distributions 



We adopt the charge density profiles from ref. pi 

Charge density profile for Oxygen sO 16 is determined by harmonic 
oscillator model 



9 \ / / \ 2 /\4 N 



where: 

p(0) = 0.141 fm~ 3 R = 1.833 fm C = 1.544. Ci = 0. 

In the case of \%Ar and the two parameters Fermi model gives 

the following profiles: 

K^ (1 + exp f_ e)/Ci)) - in, 

in the case of Argon: 

p(0) = 0.176 fm" 3 , C = 3.530 fm, Ci = 0.542 fm. 
in the case of Iron: 

p(0) = 0.163 fm -3 , C = 4.111fm, Ci = 0.558fm. 
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Ratio of effective mass to free mass for 8 16 , igAr 40 , 26 Fe'^ 
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Fi gUie 1 : The ratio of effective mass to the free nucleon mass calculated by solving of self-consistency 
equation J5j for a local value kp(r). 
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Quasielastic neutrino differential cross section on: 8 , 18 Ar 



>R Fe 



Oxygen: k F -199 MeV 
Oxygen: density profile 
Argon: k F =217 MeV 
Argon: density profile 
Iron: k F =217 MeV 
Iron: density profile 




300 400 500 600 

Transfer of energy q [MeV] 



Fi gUie 2: Comparison of differential cross sections for the scattering of neutrinos on nucleus with local 
and constant baryon density. The model of nucleus is the free Fermi gas (M* = 939MeV). Computations 
were done for neutrino energy 1 GeV. 



Quasielastic neutrino differential cross section on: 8 , 18 Ar , 2 6 F e 
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Fi gUIG 3l Comparison of differential cross sections for the scattering of neutrinos with nucleus with local 
and constant baryon density. The model of nucleus is the free Fermi gas with effective mass. Computations 
were done for neutrino energy 1 GeV. 
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Quasielastic neutrino differential cross section on: 8 16 , 18 Ar 40 , 2 6Fe ! 



Oxygen: k F = 199 MeV 
Oxygen: density profile 
Argon: k F = 217 MeV 
Argon: density profile 
rron:k F = 217MeV 
Iron: density profile 
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Figure 4: Comparison of differential cross sections for the scattering of neutrinos with nucleus with local 
and constant baryon density. The model of nucleus is by Fermi gas with RPA corrections (Af* = 939MeV). 
Computations were done for neutrino energy 1 GeV. 



Quasielastic neutrino differential cross section on: 8 , 18 Ar , 2 6Fe 



Oxygen: k F = 1 99 MeV, M = 691 MeV 
Oxygen: density profile 
Argon: k F = 21 7 MeV, M = 631 MeV 
Argon: density profile 
Iron: k F = 217 MeV, M = 635 MeV 
Iron: density profile 
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Fi gUIG 5l Comparison of differential cross sections for the scattering of neutrinos with nucleus with 
local and constant baryon density. The model of nucleus is the Fermi gas with RPA corrections and 
effective mass. Computations were done for neutrino energy 1 GeV. 
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Quasielastic neutrino total cross sections 




Figure 6; Total cross sections for scattering of neutrino on Oxygen, Argon and Iron. Calculations 
were done by applying corresponding charge densities profiles(fc.F and M* are local) and RPA. 
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